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Abstract 



We present the numerical solutions for the I — V characteristics, and 

describe the motion of fluxons in a frustrated Josephson junction made of an 

unconventional triplet superconductor and an s-wave superconductor. In the 

inline geometry, and long length limit the moving integer fluxon interacts with 

the bound fractional fluxon but is not able to change its position or polarity. 

We observe different modes of multifluxon propagation. In the small length 

limit the moving fluxon is a combination of the stable solutions that exist in 

the static case and additional steps are introduced in the I — V diagram. 
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I. INTRODUCTION 



The order parameter symmetry in the recently discovered superconductor Sr2Ru04 is 
an active area of current work. The Knight shift shows no change when passing through the 
superconducting state and this indicates that the pairing state is triplet [0. The magnetic 
field is spontaneously induced as shown by fiSR experiment, and this is a clear indication 
that the time-reversal symmetry is broken ^ . Inelastic- neutron-scattering measurements on 
single crystals of Sr2Ru04 show that the pairing state is highly anisotropic |Q. Also specific 
heat measurements support the scenario of line nodes within the gap as in high-Tc cuprate 
superconductors p. 

In Josephson junctions made of unconventional d-wave superconductors, the Josephson 
critical current can be negative depending on the orientation angles of the crystallographic 
axes with respect to the junction interface. A negative critical current can be thought of as 
a phase shift of vr at the junction interface. A half fluxon (antifiuxon) is trapped at a — vr 
junction. It's existence is confirmed experimentally by measurements of the critical current 
versus the magnetic flux in corner junctions or corner SQUID where a dip appears in the 
critical current for magnetic field equal to zero 0. The half magnetic flux quantum, (vr 
fluxon) has been directly observed using a scanning superconducting quantum interference 
device microscope in tricrystal frustrated junctions in superconducting YBa2Cu307_5 0. 

The static properties of one dimensional frustrated Josephson junctions between singlet 
and triplet superconductors have been presented for different nodal and nodeless pairing 
states 1^. The critical current and the spontaneous flux show a characteristic modulation 
with the junction orientation, which can be tested by experiment. 

The Josephson junction between s-wave superconductors supports modes of resonant 
propagation of fiuxons 0]. In the plot of the current- volt age (/ — V) characteristics these 
modes appear as near-constant voltage branches known as zero field steps (ZFS) |lTO| , p!I| . 
They occur in the absence of any external field. The ZFS appear at integer multiplies of 
Vi = hcs/2eL, where cs is the velocity of the electromagnetic waves in the junction, and L 
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is the junction length. The moving sohton is accompanied by a voltage pulse which can be 
detected at the junction's edges. 

In Josephson junctions made of unconventional ci-wave superconductors the bound half- 
fiuxon (-antifluxon) reverses its sign and emits an integer fluxon (antifluxon) when biased 
by an external current |T^-|n[. In this work we study the dynamic properties of a frustrated 
junction, between singlet and triplet superconductors and calculate the I — V characteristics. 
For the triplet superconductor Sr2Ru04 we shall assume two possible pairing states of two 
dimensional order parameter, breaking the time reversal symmetry. The first one is the 



nodeless p-wave order parameter with Eu symmetry [|T5[. The other one is the /-wave state 



proposed by Hasegawa et al, having Big x symmetry |jT6|. 

We study both the long and short length junction limit. In the long length junction 
limit, the external current cannot move the fractional fluxon (antifluxon) ff{faf) which is 
confined at a; = and the ZES exist at integer values of the dc voltage Vi. In the short 
length junction limit, additional modes in the I — V diagram are introduced at half integer 
values of the dc voltage Vi and this can be used to distinguish the possible pairing states in 
the junction. 

In the following we present the theoretical model for the corner junction in Sec. II. We 
present the results for the large junction limit in Sec. III. The case of the shorter junction 
is presented in Sec. IV and finish with the conclusions. 



II. CORNER JUNCTION MODEL 



We consider the junction shown in Fig. 1(a) between a triplet superconductor A with 
a two component order parameter and an s-wave superconductor B. The Ginzburg - Lan- 
dau (GL) equations for isotropic p-wave superconductors have been derived microscopically 
based on Gor'kov's theory [jl^ . This derivation is closely compared to the GL free energy for 
the superconducting state composed of two degenerate components. Applying the boundary 
conditions at y = and y = W where W is the width of the interface we derive the following 
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current phase relation 

J(0) = JeSin(0 + 0,), (1) 

where is defined as the relative phase difference between the two superconductors and 
0c is the intrinsic phase difference. Jc is the Josephson critical current density. For the 
type of junction that we consider where the insulator has a definite thickness and it is 
not a point contact as in the case treated by Barash et al. the Josephson effect is 
strongly directional dependent and the possibility for the tunneling of the Cooper pairs 
becomes maximum when the trajectory of the Cooper pair is vertical to the interface. So 
the total momentum of the order parameter functional depends only on the orientation 
of the interface. This dependence enters the current phase relation via the Jc-, and 0c- We 
consider a corner Josephson junction between the superconductor A and the superconductor 
B as seen in Fig. 0(b). The orientation of the a and h crystallographic axes are at right 
angles with the junction edges. We map the two segments of this junction each of length 1/2 
into a one dimensional axis. The characteristic phases 0ci, 0c2 for the two segments, for the 
pairing states that we consider can be seen in table |I[ The superconducting phase difference 
across the junction is then the solution of the sine-Gordon equation 

where 7 is the damping constant which depends on the temperature. The inline boundary 
condition reads 

where / is the inline bias current. The length x is normalized in units of the Josephson 
penetration depth given by 



A.: 



(4) 



where d is the sum of the penetration depths in two superconductors plus the thickness of 
the insulator layer. The time t is in units of the inverse of the Josephson plasma frequency 



Uq^ = \j/cs. (5) 

We can classify the different solutions obtained from Eq. ^with their magnetic flux content, 
in units of the flux quantum $o 

$ = ^(0i?-0L), (6) 

where 0l(_r) is the value of the phase at the left(right) edge of the junction. 

III. LARGE JUNCTION LIMIT 

A 4*'* order Runge Kutta method with fixed time step At = 0.01, was used for the 
integration of the equations of motion. The number of grid points is = 1000, and the 
junction length is / = 20. The damping coefficient 7 = 0.01 is used in all the calculations. 
The I — V characteristics for the first and second ZFS (corresponding to one and two 
fluxons moving into the junction) are seen in Fig. ^(a) for the Eu and in Fig. |^(b) for the 
Big X Eu pairing state. For the first ZFS two different modes of fluxon propagation exist, 
corresponding to the presence of bound fractional fluxon or antifluxon respectively at the 
junction center. For the Eu case, due to the difference in the flux content of the bound 
solutions, the // has smaller critical current than the faf. This is opposite to the Big x Eu 
case. In the Big wave case the I — V curves for the first ZFS, for the // and faf mode 



have equal critical currents ||14|| . 



For the Eu case, the external current cannot move the faf with $ = —0.75 which is 
confined at x = (see Fig. |^(a)). By applying the external current it emits an integer 
antifluxon {AF) which moves to the right and converts to an // with $ = 0.25. The AF 
hits the right boundary and transforms into a fluxon {F) which moves to the left. When the 
F reaches the center drags the // with it for a while and forms a fluxon with flux $ = 1.25 
which then breaks into a // and a F moving to the left. The fluxon hits the left boundary 
transforms into an antifluxon which moves to the center where it meets the oscillating // 
and interacts with it forming a faf and the period has been completed. A full period of 



motion back and forth takes time t = 40, and since the overall phase advance is Air, in the 
relativistic limit where u = 1 reached at high currents, the dc voltage across the junction 
will heV = 0.314. 

We may also have the situation where the // with $ = 0.25 exists at the junction center 
as seen in Fig. ^(b). By applying an external current it emits an integer fluxon which 
moves to the left and converts to a faf with $ = —0.75. The F hits the left boundary and 
transforms to an AF which moves to the right. When the AF reaches the center it meets 
the oscillating faf and forms a large antifiuxon with $ = —1.75 which then splits into a 
faf with $ = —0.75 and an integer antifiuxon moving to the right. The antifiuxon hits 
the right boundary, transforms into a fiuxon which moves to the center where it meets the 
oscillating faf, interacts with it forming a // and the period has been completed. For the 
Big X Ey_ pairing state the // has magnetic flux $ = 0.75 while the faf has magnet flux 
$ = —0.25 and an integer fiuxon or antifiuxon propagates into the junction and interacts 
with the trapped fiuxon or antifiuxon. 

The different character of the various fiuxon solutions can also be seen from the plot of the 
instantaneous voltage 0t at the center of the junction for the various fiuxon configurations. 
That plot is seen in Fig. |^ for the solutions regarding the first ZFS for the E^ state. During 
the time of one period three peaks appear in this plot by the time the fiuxon (antifiuxon) 
passes through the junction center. Note that the characteristic oscillations of (pt between 
the peaks are due to the oscillation of the bound solution about the junction center. These 
oscillations become more distinguishable in the faf case due its larger fiuxon content. Note 
also the difference in height between successive peaks in the (pt vs t diagram. This difference 
become more pronounced in the faf case due to its larger flux content. The plot of (pt at 
the edges shows two peaks during the time of one period at time instants which differ by 
half a period. So for the first ZFS the (pt vs t plot can be used to probe the existence of // 
or faf at the junction center. For the Big x Eu pairing state (not presented in the figure) 
the // has larger flux and the difference in the peak heights and the inter-peak oscillations 
become more distinguishable in the // case. 
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For the second ZFS we found four two-fluxon configurations with distinct I — V curves. 
Depending on the distance between the two vortices which is kept constant we can cat- 
egorize the solutions as seen in Fig. ^ So compared to the case of conventional s-wave 
superconductors [|1T| junction we observe several curves for the second ZFS depending on 
the relative distance between the fiuxons and this may be used to probe the presence of 
intrinsic magnetic flux. 



IV. SHORTER JUNCTION LIMIT 

We now consider the case where the junction is of relative short length / = 2. We plot 
in Fig. H the I — V characteristics for the ZFS. We see that additional steps appear at 
half integer values of the dc voltage Vi besides the ones that exist at integer values. In 
Fig. |^(a) we present (j){x), for the first ZFS, for one period (T = 4) at various instances in 
time, which are separated by At = 0.2. For the case, for t = a faf with $ = —0.75 
exists at the junction center. By applying the external current it converts to a // with 
$ = 0.25 and a faf with $ = —0.75. The faf drags the // with it forming a combination 
of fractional fluxon-antifiuxon (// — faf), with negative magnetic flux $ = —0.5 which 
moves to the right. It hits the right boundary and the fractional antifiuxon with $ = —0.75 
goes to a fractional fluxon with $ = 1.25, while the fractional fluxon with $ = 0.25 goes to 
a fractional antifiuxon with $ = —0.75, forming again a combination of fractional fiuxon- 
antifiuxon (// — faf) with total flux $ = 0.5, which moves to the left. Although the total 
flux is half integer the moving fluxon to the right direction has different structure than the 
combination that moves to the left and this explains the asymmetry between successive 
peaks in the (pt diagram seen in Fig. ^(a). The // — faf hits the left boundary transforms 
into an // — faf with $ = —0.5 which moves to the center where the period is completed. 
A full period of motion back and forth takes time t = 4, and since the overall phase advance 
is 27r, in the relativistic limit where u = 1 reached at high currents, the dc voltage across the 
junction will be = tt/2. This is indeed the value obtained from the numerical simulation 
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as seen in Fig. ^ for the solution labeled as 1/2. Note that this value is half than the case 
where a full fiuxon moves into the junction. 

In Fig. |^(c) we present (pix), for the case where a bound combination of two fractional 
vortices with total magnetic flux $ = 1.5 is propagating into the junction. In a junction of 
length / = 2 the propagating fiuxon accomplishes an overall phase advance of 6tt in a full 
period T = 4. Thus the dc voltage across the junction will be = 37r/2. The ZFS seen in 
Fig. P corresponds to a solution labeled as 3/2. Note the the structure of the fiuxon moving 
in the forward and baclcward direction is not the same and the successive peaks in the (pt vs 
t at the center of the junction do not have equal heights as seen in Fig. |^(c). 

In Fig. |^(d) we present for the case where the moving combination corresponds 

to magnetic fiux equal to 2. Thus the dc voltage across the junction will be = 27r. The 
corresponding ZFS is labeled as 2 in Fig. |^. Due to symmetry the (pt vs t at the center of 
the junction will have the same height for the forward and baclcward direction as seen in 
Fig. 1(d). 

For the Big x state for the first ZFS, by increasing the bias current, the faf with 
$ = —0.25 is transformed into a // with $ = 0.75 and a faf with $ = —1.25 forming a 
fractional fiuxon with total fiux $ = —0.5 that moves to the right. The reflected fiuxons 
has $ = 0.5 but differently to the Eu case it is composed by a // with $ = 0.25 and a faf 
with $ = —0.75. This difference can be seen in the (pt vs t diagram at the center which is 
displaced in time by half a period compared to the E^ case. However the resulting I — V 
are similar to the E^ case. 

V. CONCLUSIONS 

We analyzed the dynamics of fiuxons moving in a frustrated Josephson junction with 
triplet pairing symmetry, and calculated the I — V characteristics. The external current 
cannot move the ff{faf) which is confined at x = 0. However the external current is able 
to reorient the ff{faf) and emit an integer fiuxon(antifiuxon). For the first ZFS we found 
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two distinct curves with different critical currents, which correspond to the case where 
the moving fluxon or antifluxon interacts with a bound fractional fiuxon and antifiuxon 
respectively. The critical currents are different for the Eu and the Big x Eu pairing states 
and this can be used to distinguish the pairing symmetry. When there are more than one 
integer vortices moving in the junction, there is a possibility of different modes of fluxon 
propagation which correspond to different critical currents. 

In the small junction limit, due to the presence of the internal flux, the moving integer 
or half integer vortices will have internal structure that is formed from the combinations of 
the static solutions. The different modes in the I — V diagram exist both at integer and half 
integer values of the dc voltage. However the I — V is similar for the and Big x E^ and 
therefore it can not be used to distinguish between the two pairing states. 
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(a) 



A 




(b) 



p, or f-wave 

s-wave 

FIG. 1. (a) View of the junction between a two component triplet superconductor A and 
singlet superconductor B. The shaded region marks the interface which extends from y = 
y = W. (b) The corner junction geometry. 
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FIG. 2. I — V characteristics for the inhne geometry, for the first and second ZFS, for the 
pairing symmetry. The solutions for the first ZFS are the //, faf corresponding to a bound fluxon 
or antifluxon in the junction's center. For the second ZFS the solutions are labeled by their relative 
distance //x, where / = 20 is the junctions length and x = 1, 2, 3, 6 respectively. 
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FIG. 3. Phase 0(a;) vs x for the solutions in the first ZFS, at various instants, during one 
period separated by At = 2.8. The curves are shifted by 0.5 to avoid overlapping. Z = 20, / = 1.6, 
7 = 0.01: (a) //, (b) faf. The pairing symmetry is Eu- 
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FIG. 4. Instantaneous voltage in the middle of the junction {x = 0) vs time t, for the solutions 



in the first ZFS. Z = 20, 7 = 0.01, / = 1.6: (a) //, (b) faf. The pairing symmetry is Eu- 
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FIG. 5. Phase (f){x) vs x, for the solutions in the second ZFS, at various instants separated by 

At = 2.8. The curves are shifted by 0.5 to avoid overlapping. Z = 20, 7 = 0.01 : (a) 1/6, I = 1.6, 

(b) 1/3, I = 1.6, (c) 1/2, I = 1.6, {d)l, I = 1.6. The pairing symmetry is Ey,. 
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FIG. ^. I — V characteristics for the inUne geometry, for the pairing state, I = 2. 
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FIG. 7. Phase (p{x) vs x for the solutions, at various instants, during one period separated by 
At = 0.2. The pairing state is Eu. The curves are shifted by 0.5 to avoid overlapping. I = 2, 
7 = 0.25, 7 = 0.01. 
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FIG. 8. Instantaneous voltage in the middle of the junction {x = 0) vs time t, for the various 
solutions. The pairing state is Eu. I = 2, j = 0.01, / = 0.25. 
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